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^D ' The beliavior of Wilson loops with fractional charge is used to study the ^-dependence 






o\ 



of the free energy density e{9) for the CP^, CP^, and CP^ sigma models in two spacetime 
dimensions. The function e{9) is extracted from the area law for a Wilson loop of charge 
q = 9/2Tr. For CP^, e{6) is smooth in the region 9 ^ tt and well-described by a dilute 
instanton gas throughout the range < 6* < 27r. For CP^ and CP^ the energy exhibits a 
C^ ' clear cusp and evidence for discrete, degenerate vacua at 6' = tt, as expected from large N 

JL , arguments. For CP^ the 0-dependence is in good quantitative agreement with the leading 

^ . order large N prediction e{9) — ^Xt^^ throughout the range < < tt. 

> ■ I. INTRODUCTION 

cn 

"^ I The structure of topological charge in the QCD vacuum is central to an understanding of 

Tjj- . low-energy hadron dynamics. Lattice calculations have provided quantitative confirmation of the 

o ■ 

00 ' Witten-Veneziano [1|] relation between the topological susceptibility of quenched QCD and the mass 

o ■ n 

of the rj' meson [4]. The success of the Witten-Veneziano relation, which relates the topological 
k>( \ susceptibility to the chiral structure of the quenched rj' correlator |3|] is most easily understood in 

j^ ■ the context of the large A^ approximation. (In this case, the large N approximation provides a 

justification for resumming quenched "hairpin" insertions.) The observation of extended, coherent, 
codimension 1 sheets of topological charge in Monte Carlo generated SU(3) gauge configurations 
[^,|5| points to a new paradigm for the QCD vacuum. Contrary to the standard instanton liquid pic- 
ture, the vacuum is apparently a "topological sandwich" consisting of alternating sign membranes 
of topological charge. Longstanding arguments of Witten [g] based on large-N chiral dynamics 
show that, for sufficiently large N, instantons should disappear from the QCD vacuum and be 
replaced by codimension-1 membranes which are in fact domain walls between discrete "k-vacua," 
where the effective local value of 9 differs from that of the true ground state by Ink, with k an 
integer. The topological sandwich picture which has emerged from the Monte Carlo calculations 



is quite compatible with these large N arguments. This picture is in some ways a 4-dimensional 
analog of Coleman's picture of ^-dependence in the massive Schwinger model [3]. In that model 
there are no instantons, and 6 can be interpreted as a background electric field. The domain 
walls are just charged particles, whose world lines have codimension 1. These world lines separate 
regions in which the electric flux differs by one unit, i.e. A6 = ±2tt. A similar interpretation of 
^-dependence applies to the two-dimensional CP sigma models, which also have a U(l) gauge 
invariance. Luscher J8«] clarified the analogy between 2D U{1) theories and 4D Yang-Mills theory 
by pointing out the similar role played by the Chern-Simons tensors in the 2D and 4D theories. 
The Wilson loop in the CP model is analogous to a "Wilson bag" in QCD, which is an integral 
of the CS tensor of the Yang-Mills field over a three-dimensional surface. From this viewpoint 
the Wilson loop <f A ■ dl in 2D should be interpreted as a one-dimensional surface integral of the 
Chern-Simons flux e^^yA'^ . Monte Carlo studies of the CP models [9|, |lO(] have shown that for 
A^ > 3 these models exhibit a topological structure quite analogous to that observed in lattice 
QCD, with the vacuum occupied by extended coherent topological charge structures of codimen- 
sion 1 [y]. In striking contrast, the topological charge distributions in the CP^ and CP^ models 
are found to be dominated by small instantons [10]. The CP^~^ models thus provide a detailed 

example of Witten's arguments that instantons should "melt" or become irrelevant at large N . For 

I I 

the CP models, the instanton melting point is found to be about iV ~ 4 [10]. (As discussed in 



Hi, the value of A^ at which instantons become irrelevant may depend on the lattice formulation, 
but for any latticization, at least the CP^ model is expected to be dominated by small instantons.) 
Although the value of the 9 parameter in real-world QCD is zero to high accuracy, a deeper 
understanding of topological charge structure can by obtained by considering gauge theories with 
nonzero 6 terms. Unfortunately, Monte Carlo studies of QCD at nonzero 9 are severely hampered 
by the fact that the 6 term contributes an imaginary part to the Eucidean action. As a result, the 
exponentiated action in the path integral cannot be interpreted as a probablility, precluding the 
direct simulation of such theories by Monte Carlo techniques. Naive reweighting methods, in which 
the e'^'' factor in the path integral is introduced in the ensemble average over ^ = configurations, 
are sufficient for probing small values of 9. But some of the most interesting issues associated with 
multi-phase structure are most directly addressed by studying the region ~ vr, where reweighting 

techniques for extending the reach of Monte Carlo 



12 



13, h4, UM, but definitive results near 9 ^ ir are 



methods are ineffective. More sophisticatec 
techniques to large 9 have been introduced 
still difficult to obtain by these methods. 

In this paper we present a method for exploring ^-dependence in two-dimensional U{1) gauge 



theories which is based on the calculation of fractionally charged Wilson loops in the 6 = theory. 
The results for the free energy density e{6) = E{6)/V in the CP models exhibit the power 
of this method, providing direct numerical evidence for a first order phase transition at 9 = ir for 
the CP^ and CP^ models. In clear contrast with the larger A^ models, the instanton dominated 
CP^ model exhibits smooth behavior in the 9 = it region, as expected from a dilute instanton gas 
calculation. (Our results do not rule out a second order phase transition at = tt in the CP^ 

n 

model, which is expected from theoretical arguments [16]). 

II. FRACTIONALLY CHARGED WILSON LOOPS 

The central observation which we exploit here is that, in the path integral for a 2D U(l) gauge 
theory, including a closed Wilson loop with charge q = 9/2ir is equivalent to including a 9 term 
in the two-volume enclosed by the loop. The Wilson loop goes around the boundary dV of a 
two- volume V, 

9 [ (fxQ{x) = — i A-dx (1) 

Jv 27r Jqv 

where 

Q{x) = ^e.^F^^" (2) 

is the topological charge density. Since the CP^^^ models have nonzero topological susceptibility, 
the free energy per unit volume £{9) inside the loop is greater than e(0) outside. This gives rise to 
a linear confining potential between fractional charges or equivalently, an area law for large Wilson 
loops. The coefficient of the area law determines the difference in vacuum energy inside and outside 
the loop: 

{Wc{q)) ~ exp [- {e{9) - ^(O)) V\ \e=2., (3) 

where C = dV . Note that since e{9 + 2tt) = £{9), the coefficient of the area law is periodic in the 
charge q. At integer values of charge, the confining force is completely screened and the area term 
in the Wilson loop vanishes. The behavior observed in the region q ~ 1/2 distinguishes between an 
instanton-dominated model and a large A^ domain wall model. For an instanton theory, a dilute 
gas approximation gives 

e{9)-e{0) = xt{l-cos9) (4) 




FIG. 1: Behavior of the topological susceptibility as a function of N, divided by the large N prediction, for 
both log plaquette and overlap definitions of lattice topological charge. 



where xt is the topological susceptibility. On the other hand, large N considerations [Gj] predict 
quadratic behavior, with periodicity leading to cusps at odd multiples of vr. 



e{e) - e(0) = ^xt miukez {0 - 27rkf 



(5) 



Physically, these cusps represent "string-breaking" or vacuum screening, which occurs when it is 
energetically favorable to screen the background electric flux by one unit from ^ + 6to—2+S. For 
large loops, this can be interpreted as a first order phase transition taking place inside the loop. 

We study the Wilson loops for the CP^, CP^, and CP^ models. As shown in Ref. [lol, CP^ 
is instanton-dominated, while CP^ and CP^ are both above the instanton "melting point," It 
was found that the topological charge distribution for CP^ is dominated, at large /3, by small 
instantons, while for CP^ and CP^ instantons do not occur but topological charge instead appears 



in the form of coherent one-dimensional domain wall- like structures. In Ref. 
topological suscepibility were compared with the leading-order large N prediction 

3^2 



lOJ I results for the 



Xt 



AirN 



(6) 



where ^ is the mass gap (nonsinglet meson mass). As shown in Fig. 1, the result for CP^ is still 
significantly above the large A^ prediction, but CP^ is in good agreement. Fig. 1 plots the ratio of 
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FIG. 2: The free energy density e{9) for CP^ (x's), CP^ (D's), and CP^ (O's) as measured from fractionally 
charged Wilson loops. The lower and upper curves are the instanton gas and large N predictions, normalized 
to the same topological susceptibility. Note that, for 9/2tt > ^ the large N curve represents the energy of 
the false (unscreened) vacuum. 



the measured value of xt to the large A^ prediction ([5]). The x's and o's denote the results for log 
plaquette and overlap topological charge, respectively. For each of the CP models, the value 
of the coupling constant /? was adjusted to give a correlation length of « 5 to 7. [The CP^, CP^, 
and CP^ data were obtained at /? = 1.2,0.9, and 0.8 respectively, for which the meson masses are 
fi = 0.179(3), 0.186(3), and 0.212(2). Note that for CP^ and CP^, xt scales properly as xt/i^ while 
for CP^ the small instantons cause Xtl l^^ to diverge in the continuum limit.] The calculations were 
done on both 50 x 50 and 100 x 100 lattices and the effect of finite volume was found to be negligible. 
The results plotted in Figure 1 indicate that by studying these three models, we cover the entire 
range of topological charge dynamics from CP^ which is instanton-dominated, to CP^ where bulk 
topological properties are described with reasonable accuracy by the large N approximation. This 
assertion is further supported by the results for dependence of the vacuum energy presented in 
this paper. 

One of the most striking features of the Monte Carlo results for fractionally charged Wilson 
loops is the difference between the behavior in CP^ versus that in CP^ and CP'^ in the region 



IT < 9 < 2tt. Since £{9) is periodic and an even function of 0, its value in this region should be 
determined by its value in the range < < vr by reflection around vr, 

eie) = e(2^ - 9) (7) 

As seen in Figure 2, the measured value of e{9) for CP^ (^'s) is, within errors, nicely periodic and 
symmetric around 9 = ir, and in fact fits well to the dilute instanton gas prediction (j3|) throughout 
the range < 9 < 2ir. On the other hand, for CP^ (D's) and CP'^ (o's), as shown in Fig. 2, the 
coefficient extracted from a simple area-law fit to the Wilson loops continues to rise beyond 9 = n, 
violating the expected symmetry ([7]) . We will argue that the behavior of CP^ and CP^ for ^ > vr is 
an effect of having two nearly degenerate ground states. This behavior can be easily understood in 
terms of the large N picture[6[], in which there are two nearly degenerate quasi- vacua when ^ ~ vr. 
One vacuum has a background electric field 9/27: ~ + 2 • The other is the one in which a unit of flux 
has been screened, so that 9/27r ~ — 2- A Wilson loop with length R in the spatial direction and 
T in the time direction can be interpreted as the T-dependent propagator of a "string" of length 
R, consisting of a +q and a —q charge with an amount q = 9/271 of electric fiux between them. 
This state has a large overlap with the vacuum state containing background fiux of 9 /27x. But for 
9 > 7T, the true ground state is the one where the flux has been screened by one unit to 9/27T — 1. In 
order for the Wilson line to couple to this screened vacuum, the fiux string must break via vacuum 
polarization. It is thus expected that, for 9 > 7r, the Wilson line will have a much larger overlap 
with the false (unscreened) vacuum than with the true (screened) vacuum. As a result, the Wilson 
loop area law tends to be determined by the energy of the unscreened vacuum, even for 9 > 77 
where the screened vacuum has lower energy. Since the Wilson line couples preferentially to the 
unscreened vacuum, we expect that our results for £{9) are measuring the true ground state energy 
throughout the range < ^ < vr, where the unscreened vacuum is the true vacuum. By invoking 
the refiection symmetry ([7]) we obtain a complete determination of £{9). 

The results for e{9) for CP^ are plotted for < < vr in Fig. 3. With the topological sus- 
ceptibility fixed to the value obtained from the fiuctuation of the integer- valued global topological 
charge, xt = (^^)/^i the solid line is a zero-parameter fit to the dilute instanton gas formula ^. 
Also plotted (dotted line) is the leading order large N prediction, 

e(9)-eiO) = ^Xte' (8) 

The corresponding results for CP^ and CP^ are shown in Figs. 4 and 5. The solid and dotted 
curves are again the instanton gas and large A^ predictions with normalization fixed to the measured 
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FIG. 3: e{e) for the CP'^ model. 



■i? 3 
LU 



o 



0.1 



0.2 0.3 

e/2jt 



1 1 1 


1 1 


- 


,• ' 1 






;c 










; 






,' 










" 






< 


< 




; 




: 


< 












F'' .^^-^^'^'^^^^^ 


.■' > 


^^^^ 


— — —""''''''"^ 1 1 1 


1 1 



0.4 



0.5 



FIG. 4: s{9) for the CP'^ model. 




FIG. 5: e{e) for the CP^ model. 

Xt- Unlike the CP^ case, the data for CP^ and CP^ shows a clear departure from the instanton 
gas model in the direction of large N. In particular, there is convincing evidence that e{9) has a 
positive slope at ^ = vr, and since e{9) = e{2TT — 6), this implies a cusp with discontinuous derivative 
at 9 = n. Moreover, the results for e(^) for the CP^ model (Fig. 5) are in good agreement with 
the large N prediction ([5]) in both the magnitude of xt and in the shape of the function e{9). 
(The Xt used here is from the log plaquette charge, but for CP^ there is little difference between 
the log plaquette and overlap values for xt- See Fig. 1.) In all of the models, the value of xt 
obtained from fractionally charged Wilson loops using xt = ^"(O) is in excellent agreement with 
that obtained from the fluctuations of the global topological charge xt = (i^^)/^- This agreement 
provides additional evidence that the fractionally charged Wilson loop method gives an accurate 
determination of e{9) over the entire range < ^ < vr. 

III. LATTICE CP'^-^ MODELS 



The Lagrangian for CP sigma models in the continuum is 



L = d^z*d^Zi + {z*d^Zi){z*d^Zj) 



(9) 



where z* are N complex fields, i = 1,. . . ,N, satisfying a constraint z*Zi = 1. This Lagrangian is 
invariant under a local U(l) gauge transformation: Zi[x) — > e^^^^' Zi{x), for arbitrary space-time 
dependent a{x). We can introduce a dummy gauge field A^ and rewrite the Lagrangian as 

L = {d^-iA^)z*{d,, + iA^)z^ (10) 

where 



An = -i 



zlid,z)-id^zl)z,] (11) 



To put CP models on the lattice, we introduce U{1) link fields U{x,x + fi) = e f*-^^. 
CP fields are defined on the sites as Zi{x). The lattice action consists of gauge- invariant 
nearest-neighbor hopping terms, 

S = l3N^Zi{x)'^U{x, X + fi)zi{x + /t) + c.c. (12) 

This lattice action is used in our Monte Carlo simulation to generate an ensemble of field configu- 
rations. The Zi fields are updated by a Cabibbo-Marinari heat bath algorithm (i.e. by applying an 
SU{2) heat bath to all pairs of components Zi, Zj), while U{1) link fields are updated by a multi-hit 
Metropolis algorithm. 

In the continuum theory, the topological charge density is proportional to the electric field 
strength, 

Qi^) = ^e^ud'^A'^ (13) 

Zn 

The simplest definition of the local topological charge density on the lattice is thus obtained from 
the plaquette phase, 

Q, = ^lnP, (14) 

zvr 

n recent studies of topological charge structure in the CP^~^ models[9|, [lOj] (as well as in QCD 
J]), it was found that a "fermionic" construction of the lattice topological charge density based on 
the overlap Dirac operator provides a much clearer view of coherent structures (both domain walls 
and instantons). The density obtained from the plaquette phase has a high level of short range 
anticorrelated noise which is smoothed out in the overlap-based topological charge distribution. 
Nevertheless, the two choices of lattice topological charge give similar results for bulk topological 
susceptibility, particularly for larger N (see Fig. 1). 
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FIG. 6: Static potential for the CP^ model with Wilson loop charge q — 0.3. 



The Wilson loops we consider in this paper can be interpreted as world lines of fractionally 
charged test particles. In the compact formulation we are using, the link phase associated with a 
particle of charge g is V^, a fractional power of the link variable. In general this would require a 
branch choice which will not be gauge invariant. However, for U{1) gauge theories in two space- 
time dimensions, there is a straightforward gauge invariant procedure for calculating closed Wilson 
loops with fractional charge, using the lattice Stoke's theorem, 

(W(C)) = (HC^/) =exp(z^lnP, (15) 

leC ueV 

where the loop C is the boundary of the enclosed volume V, I labels the links on the loop and 
u labels the plaquettes enclosed. Here P^ is the gauge invariant product of links around a single 
plaquette, with its log taken to be on the principal branch, 



vr < — i In P,, < vr 



(16) 



(Note that this branch choice is gauge independent, unlike a branch choice for the link phase V^.) 
To extract the ground state energy density e{6) from the Wilson loop expectation values, we 
consider rectangular loops Wq{R,T) with sides R and T and charge q. The fitting of the Wilson 
loops for charges q < ^ is particularly simple in that there is very little deviation from a pure area 
law, even for loops as small as 2 x 2. This allows a surprisingly accurate determination of e{6) from 
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FIG. 7: Static potential for the CP^ model with Wilson loop charge q — 0.3. 
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FIG. 8: Static potential for the CP^ model with Wilson loop charge q — 0.3. 
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the area law coefficient over tfie entire range < ^ < vr. It is convenient to extract the area law 
coefficient by first studying the static potential Vq{R). Regarding T as Euclidean time, the static 
potential is calculated by fitting the data at each fixed value of R to an exponential T dependence, 

{Wq{R,T)) ~ const. X exp(-yg(i?)r) (17) 

The T-dependence for any fixed R is well described by a single exponential for T > 1. Moreover, 
the static potential Vg{R) is found to be quite accurately linear, 

Vg{R) ~ i? X e{2TTq) (18) 

The static potential for q = 0.3 is shown in Figs. 6, 7, and 8 for CP^, CP^, and CP^, respectively. 
The value of e{9) is determined by fitting the slope of the linear potential. 

As we have discussed, for the CP^ model, the function e{9) decreases for 6 > tt and is nicely 
consistent with the expected symmetry ([7]). At 6 = 27r, the slope of the linear potential for CP^ 
is consistent with zero, as seen in Fig. 9. For the CP^ and CP^ models, the behavior for 6 > tt 
is quite different. As we have discussed, the measured value of e(^) continues to increase beyond 
6 = TT and violates refiection symmetry around vr. The explanation of this behavior in terms of 
coupling to the unscreened false vacuum suggest that, for g > i, we might be able to see evidence 
of the true vacuum by observing a levelling off of the slope of the linear potential for very large 
values of R (and T) as the effect of vacuum screening sets in. Within the limits of our statistics, we 
have not been able to see significant evidence for a change of slope in V{R), even for loops as large 
as 10 X 10. Figs. 9 and 10 show the results for V{R) with a Wilson loop charge of 1.0. We see that 
the CP^ and CP^ results in Fig. 10 still indicate a linear confining potential. The situation here 
is analogous to the problem of observing string breaking in QCD, i.e. the screening of the linear 
potential by closed quark loops. The direct observation of string breaking from Wilson loops in 
full QCD is notoriously difficult (c.f. [ilT,]), and most studies of this process have had to resort to 
coupled channel calculations. 

IV. INSTANTONS VS. LARGE N 

The 6 dependence of the free energy that we have calculated from fractionally charged Wilson 
loops provides a detailed view of the transition between instanton physics and large-N physics in 
an asymptotically free gauge theory. The ability to obtain accurate results for the free energy 
as a function of 9 over the whole range < 9 < tt opens new possibilities for the Monte Carlo 
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FIG. 9: Static potential for the CP^ model with Wilson loop charge q = 1.0. The solid line is the best 
zero-slope fit. 




FIG. 10: Static potential for the CP^ and CP^ models with Wilson loop charge q = 1.0. 
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exploration of topological phase structure in two-dimensional gauge theories. The understanding 
of this phase structure has broad implications, not only as an analog of QCD vacuum structure 
in four dimensions, but also for the structure of topologically ordered states in two-dimensional 
superconductors [18] . The appearance of a cusp in e{9) at 9 = tt indicates the presence of a first 
order phase transition, pointing to the existence of two discrete vacuum states which are crossing 
in energy at ^ = tt, as first suggested by Witten [6] for both QCD and CP models. The long 
range coherent codimension 1 topological charge membranes which have been observed in Monte 
Carlo configurations \4, Isl are naturally interpreted as domain walls between discrete vacua. 

It is interesting to consider the calculation of e{9) in the dilute instanton gas approximation 
@ and compare it with the assumptions that lead to the large N prediction ([5]). We write the 
Euclidean path integral with periodic boundary conditions as a sum over "winding number" (global 
topological charge) sectors, 



Z= Yl ^•^^"'^ (^9) 



where Zi, is the path integral over configurations with total topological charge i^. For a dilute 
instanton gas, this is given as a sum over instantons and antiinstantons. 



z^i z"^ 



Zu= Yl —,^5u,u,-u,=Iu{2z) (20) 

Here z is the partition function for a single instanton or antiinstanton in a box of volume V. Thus 
the partition function is 

oo 

Z= Y Iu{2z)e'''^ =exp{2zcos9) (21) 

This gives the result @ for the free energy e{9) = — InZ/V. 

To contrast this result with the large N prediction ([5]) we go back to the sum over topological 
charge sectors and employ a Poisson resummation, 

/oo ^ 

dQZQe'Q'Y.^(Q-^^= Y. dQZge''^'^'-^^''^ (22) 

Here we assume that for large volume, the partition function over topological charge sectors Z^ can 
be smoothly interpolated between integer values Q = i^ by a continuous function Zq. We can now 
show that the 9 dependence ([5]) predicted by large N arguments actually follows more generally 
from the assumption that the fluctuations of global topological charge are gaussian. In particular, 
the integer k that labels discrete k- vacua is dual to the topological charge v in the sense of Poisson 



Z = const. X ^ expl ^{e-2TTkf\ (24) 
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resummation. If we take 

Zq = Zoe-«'/2>^*^ (23) 

then the partition function can be written 

Vx 
exp ^ - 

A;=— oo 

Note that the Poisson resummation ()22p has taken us from a sum over global topological charge 
sectors, in which each term is periodic in 9, to a sum over discrete /c-vacua, where periodicity in 9 
is obtained by the fact that the k th term becomes the (/c + 1) th term when 9 — > 9 + 2tt. Upon 
taking the logarithm to get the free energy, a single term in the sum (I24p . with the minimum 
value of {9 — 2irk)'^ dominates the logarithm in the limit V — s- oo. This reproduces the large N 
result dS]). We see that the purely quadratic dependence of e{9) which is obtained in the large N 
approximation follows simply from the assumption of gaussian fluctuations of global topological 
charge. Conversely, the deviation from quadratic ^-dependence exhibited by the instanton gas 
model can be related to the deviation of (j20p from a pure gaussian distribution. The fact that a 



pure gaussian distribution in i^ implies a first order transition at = tt has been discussed in |l2l |. 



V. CONCLUSIONS 

The results presented here demonstrate the utility of fractionally charged Wilson loops as a 
probe of topological structure in f/(l) gauge theories. Applied to the 2D CP models, this 
method clarifies the nature of the transition from an instanton-dominated vacuum at small N to 
a domain-wall-dominated vacuum at large A^. Although in two dimensions both types of vacuum 
lead to confinement of fractional U{1) charge, the confinement mechanisms are quite different in 
the two cases. This distinction is of interest for investigating the relationship between topological 
charge and confinement in higher dimensional theories like QCD. In two dimensions, instantons 
give rise to electric charge confinement by introducing a phase incoherence in the fractionally 
charged loop arising from the varying number of instantons minus antiinstantons inside the loop, 
producing an area law falloff of the Wilson loop. Because the instanton topological charge is locally 
quantized, this phase incoherence automatically disappears for an integer charged loop. This is 
why the e{6) for the CP^ model returns to zero at 9 = 2tt and exhibits no "false vacuum" effect. 
Because topological charge is not locally quantized in the large A^, domain wall dominated models, 
the amount of charge inside the loop is not restricted to be an integer. This gives rise to phase 
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incoherence and an apparently nonvanishing area law coefficient, even for integer charged loops. As 
is well-known, the instanton confinement mechanism does not extend to higher dimensions. If quark 
confinement in 4D QCD is to be associated with topological charge fluctuations, the instanton gas 
or liquid model does not provide such a connection. On the other hand, the domain wall structure 
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and observationally \4, ly, 19|, [10|] quite 



y,y,u 



of the large A^ CP vacuum is both theoretically 
analogous to that of 4D QCD. The disordering of Wilson loop phases by layered, codimension 1 
domain walls is likely to be similar in the two cases. The study of discrete vacua and domain 
walls in CP^"^ models could possibly provide insight into quark conflnement and topological 
phase structure in QCD. For the 2D CP models, fractionally charged Wilson loops provide a 
promising tool for studying this structure. 

This work was supported in part by the Department of Energy under grant DE-FG02- 
97ER41027. 
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